Fluctuations and stochastic processes in one-dimensional many-body quantum systems 
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We study the fluctuation properties of a one-dimensional many-body quantum system composed 
of interacting bosons, and investigate the regimes where quantum noise or, respectively, thermal 
excitations are dominant. For the latter we develop a semiclassical description of the fluctuation 
properties based on the Ornstein-Uhlenbeck stochastic process. As an illustration, we analyze the 
phase correlation functions and the full statistical distributions of the interference between two one- 
dimensional systems, either independent or tunnel-coupled and compare with the Luttinger-liquid 
theory. 
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Measurement of fluctuations and their correlations 
yields important information on regimes and phases of 
many-body quantum systems In ultracold atomic 
systems, these correlations revealed the Mott insulator 
phase of bosonic Q and fermionic Q atoms in opti- 
cal lattices, they allowed detection of correlated atom 
pairs in spontaneous four-wave mixing of two colliding 
Bose-Einstein condensates [3| and Hanbury-Brown-Twiss 
correlation for non-degenerate metastable 3 He and 4 He 

]. Furthermore, they have 



and have been employed 



atoms Q and in atom lasers 
allowed studies of dephasing 
as noise thermometer [1, [§| . 

A key question in the physics of many-body quantum 
systems at finite temperature is how much of the ob- 
served fluctuations and their correlations are fundamen- 
tally quantum, and which are caused by the thermal ex- 
citations in the system. In the present Letter we address 
this problem starting from a description of the excita- 
tions in the system and their occupation numbers. This 
allows us to directly explore the contributions of quan- 
tum (ground state) noise and thermal excitations. 

We consider a quantum degenerate spin-polarized gas 
of bosonic atoms in an extremely anisotropic trap, with 
transversal confinement frequency uj± much larger then 
the longitudinal confinement frequency uiu (typically 
> 1000). If both the temperature T and the 
mean-field interaction energy per atom are small com- 
pared to the radial confinement (k^T <C Tuo±, niiya s <C 
1, where nro is the linear atom- number density and a s is 
the atomic s-wave scattering length), the atomic motion 
is confined to the radial ground state of the trapping po- 
tential. In this ID regime a "quasi-condensate" emerges, 
which can be characterized by a macroscopic wave func- 
tion with a fluctuating phase Io|- 12|. 

The statistical properties of the fluctuating phase are 
the focus of this study. They can be probed by interfering 
two identically prepared ID systems by creating quasi- 
condensates in two parallel, identical traps [HI]. When 
released they expand freely, overlap and interfere. The 



local phase of the interference reflects the fluctuating rel- 
ative phase 9{z) of the quasi-condensates. The fluctua- 
tions in 9(z) manifest themself in the phase correlation 
function C e (z - z') = (exp[id(z) - i0(z')}). The full dis- 
tribution function of the interference contrast has been 

Hill. 



derived in Refs. 



In our investigation we consider the general case of 
two ID quasi-condensates which can be tunnel-coupled 
to each other, described by the effective Hamiltonian [l6| 



: ]+ l^(z)tf(z) 



H = f dz{^2[i>}(z)T^( 



(i) 



Here HJ is the tunnel-coupling matrix element, T = 
— [ft 2 /(2m)]d 2 /dz 2 — fl with the chemical potential fl = 
gniD — M, and g = 2Tiuj^a s is the atomic ID interaction 
strength in the limit nrDa s <C 1 17 1. 

We study this system based on the description 
of the quasi-condensate properties by a spectrum of 
Bogoliubov-type modes |l8|, which are free-particle-like 
in the short-wavelength limit and phonon-like in the long- 
wavelength limit [19(. We model the "experimental" re- 
alizations of the atomic quasi-condensate fields by im- 
plementing a numerical scheme for generating the ini- 
tial conditions in the truncated Wigner approximation 



2Ct l2lj and represent their wave functions as tpj (z) 



bj and Snj 



-y/^iD + 6n,j(z) exp[i(j)j(z)], j = 1,2, where 
are the local phase and density fluctuations, respec- 
tively. We decompose these fluctuations into waves corre- 
sponding to elementary excitations of two coupled quasi- 
condensates by means of an extension of the approach 
by Mora-Castin [l8[ as developed by Whitlock and Bou- 
choulc [16]. The amplitudes and node positions of these 
waves are chosen randomly, assuming the Bose-Einstein 
statistics of the elementary excitations. In particular, the 



2 



relative phase 8(z) = 4>i{z) — 4>2(z) is modeled as 



,1/4 



#(2) = \Z2/(ni D £ max ) [(rjk + 2 gnxx>)/r)k\ ' v 



v /0 fe |ln&|sin(fcz + 27r&) 



(2) 



where are random numbers (obtained by a pseudo- 

random number generator) uniformly distributed be- 
tween and 1, i]k = {h~k) 2 /(2m) + 2hJ and the sum- 
mation is taken over the discrete spectrum of wave vec- 
tors k equal to (both positive and negative) multiples of 
27r/L max [jgj . A similar expansion holds for the density 
fluctuations. The explicit dependence of the wave am- 
plitude on the random numbers reflects the statistics 
of the occupation numbers of the elementary-excitation 
modes. To include both thermal and quantum fluctua- 
tions (zero-point oscillations of the atomic field) : 



B fe = 2- 1 coth[ £fc /(2fc B T)], 



(3) 



where £k = ^/ilkivk + 2 gnyo) is the energy of the ele- 
mentary excitation. We refer the use of Eq. (|3]) to as the 
"full Bogoliubov approach" . Alternatively, if we choose 
to neglect the quantum fluctuations: 



B k « k B T/e k . 



(4) 



We first analyze the full Bogoliubov approach and cal- 
culate the full distribution function of the contrast fol- 
lowing [li], EH : The interference pattern integrated over 
a sampling length L is characterized by the complex am- 
plitude operator A(L) (22j |. Each experimental run yields 
probabilistically a complex value A(L). The expectation 
value (A(L)) is zero, but (i(L)+i(L)) = (\A{L)\ 2 ) is not. 
It is convenient to study the statistical distribution W(a) 
where a = \ A(L)\ 2 / (\A(L)\ 2 ) is the square of the absolute 
value of the integrated contrast scaled to its mean. 

Using Eq. ((3J we generated integrated contrast dis- 
tributions for a wide range of parameters. In Fig. 1 we 
display W(a) as a function of the sampling length, both 
for zero and for non-zero coupling. The calculations have 
been done for L max = 100 /im and 60 modes taken into 
account (increasing of the number of modes to 120 does 
not change the results significantly). 

In the special case of zero tunnel coupling between the 
condensates (J = 0), statistical independence of fluctua- 
tions in each quasi-condensate allows to separate corre- 
lations: 



dz / dz' {^l(z)Mz')) (^l{z')^ 2 {z)), 



(\A(L)\ 2 ) 



and a general formula for the computation of all the mo- 
ments of W(a) can be found from the Luttinger-liquid 
formalism [lj, EEI- The stochastic properties of W(a) 
are then determined by a single dimensionless parameter 
ktL, where kt = mk^T / \h 2 hid) is the inverse thermal 
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FIG. 1: (Color online). Interference contrast distribution 
W(a) of a 87 Rb quasi condensate with niD = 59 /im and 
u)± — 2ir x 3 kHz as a function of the sampling length L 
for (a) T = 31 nK, J = 0, (b) T = 60 nK, J = 0, (c) 
T = 60 nK, J = 2tt x 1 Hz, and (d) T = 60 nK, J = 2tt x 3 Hz. 
Lines represent the results of the full Bogoliubov modeling for 
L = 10 fxrn (red), 24 fim (green), 37 /im (blue), and 51 /zm 
(black). Symbols (circles, squares, diamonds, up triangles) of 
the respective colors show the results of the Luttinger liquid 
approach 9J for the same values of L and J = 0. 



coherence length, m is the mass of the atom. For 87 Rb 
k t w 1.815 j tim- 1 (T/lQ0 nK)(10 nm^/nm). There is 
very good agreement between the full Bogoliubov cal- 
culations and the Luttinger liquid formalism, and one 
observes (Fig. la,b) the characteristic change between a 
Gumbel-like distribution to an exponential distribution 
as the ratio of the averaging length to the characteristic 
phase- coherence length grows Q 

If J ^ and large enough (i.e. J ~ 27r x 1 Hz for 
the typical experimental range of niD, T, and L), we 
observe a different picture: the distribution W(a) stabi- 
lizes at some peaked shape and preserves this shape as 
L grows further, (Fig. lc,d). This is characteristic for 
the phase locking between the two matter waves. Since 
the Luttinger liquid approach [l4|, 15 1 is based on the 



assumption of statistical independence of fluctuations in 
the two quasi-condensates, it can not easily be extended 
to the tunnel-coupled systems described by Eq. ([T]). 

We can now study the effect of quantum fluctuation by 
using Eq. (|4| instead of Eq. ([3]) . For weakly interacting 
ID systems the differences are small (Fig. 2). 

This observation suggests a simple semiclassical de- 
scription of the of the noise properties at distances longer 
than the healing length £h = ^C/(7t71id) (Luttinger pa- 
rameter JC — nhy niu / (mg)) , where density fluctuations 
are suppressed and the main contribution to noise comes 
from fluctuations of the relative phase 9(z). For ther- 
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FIG. 2: (Color online). Interference contrast distribution 
W(pt) for (a) ktIj — > oo and (b) ktIj = 1.0. Solid lines 
display the results of the Ornstein-Uhlenbeck stochastic pro- 
cess modelling for htL = 1.85 (red) and 4.43 (black). Results 
of the full Bogoliubov [Eq. |(3|] simulations and Bogoliubov 
simulations without quantum noise [Eq. @] are shown by 
triangles and circles, respectively. 



mal excitations, the fluctuations of ( 
their autocorrelation function is lit 



(z) are Gaussian and 



{6(z)8(z')) = KT ljexp(-\z - z'\/lj), 



(5) 



with I j = \^h/{mJ) . For 87 Rb lj « 5.367 ^m y/j^JJ, 
where J\ = 27r x 1 Hz. Eq. ([5]) is valid under the follow- 
ing assumptions: (i) we can neglect atom shot noise [23j . 
(ii) the density fluctuations are suppressed, (iii) quan- 
tum fluctuations of the phase can be neglected, and the 
mean occupation number for the given mode is taken in 
the classical (Boltzmannian) limit Eq. ((4]). The relative 
phase evolution along z can be described by an Ornstein- 
Uhlenbeck stochastic process 24j, where the coordinate 
z plays the role of time: 



±B(z) = -!*(*) + /(*). 



(0) 



the 



properties 
l 



Here f(z) is the random force with 
(f( z )) = 0, (f{zi)f{z 2 )) = 2k t 6(z 1 - z 2 ), and lj L plays 
the role of the friction coefficient. The local variance of 
the relative phase should not depend on z and the initial 



value 



is distributed according to a Gaussian with 



zero mean and variance (9 2 (z')) = KtIj, i.e., the station- 
ary distribution following from Eq. (|6|). 

This leads to a very simple and efficient way to calcu- 
late the fluctuation properties. We propagate 9(z) from 
z = to z = L using an exact updating formula for Eq . 
(j6]) |25j and compute for each run the complex phase [26j . 
A statistical analysis of the full distribution function of 
a on the ensemble of runs shows a very good agreement 
between the Bogoliubov simulations and the stochastic 
process modeling (Fig. [2]). 

For which parameters and observables are the funda- 
mental quantum fluctuations in ID systems observable? 

We first analyze the modification of the full distribu- 
tion function W(a). The contribution of quantum noise 
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FIG. 3: (Color online), (a) Interference contrast distribution 
W{a) for 87 Rb atoms, did = 60 /im _1 , J = 0, and L = 10 /jm 
calculated for both thermal and quantum fluctuations (solid 
line) and thermal fluctuations only (dashed line) for two dif- 
ferent parameter sets: T = 10 nK, uj± — 2tt x 3 kHz (I, red) 
and T = 40 nK, uj± = 2iv x 60 kHz (II, blue). The inset 
shows the temperature dependence of \C\ 2 = (\A(L)\ 2 ) / L 2 
(the average square of the dimensionless coherence factor) for 
the respective cases, (b) Ratio of the centered sth order mo- 
menta R s for s — 2 (solid line), s = 3 (dashed line; hardly 
distinguishable from s — 2 in the case I) and s = 4 (dot- 
dashed line) as a function of temperature for two parameter 
sets at L — 10 /im. Quantum fluctuations manifest them- 
selves in R s < 1. (c) The same as in (b) but for L = 20 fim. 
The effects of atomic shot noise on {\A(L)\ 2 ) and W(a) arc 
not taken into account. 



will be detectable in W(a) at very low temperatures 
k-&T <c /j, and short length scales L « 10 /im (Fig. 3). It 
can be quantified by the ratio R s = ((a — l) s )th/((a — 
l) s ) q where the averages (. . .) q and (. . .)th are obtained 
by either the full Bogoliubov approach including quan- 
tum fluctuations [Eq. ([3"])] or considering only thermal 
fluctuations [Eq. respectively (Fig. 3b, c). If the 

quantum noise is negligible, then R s = 1. 

We can quantify the relative contribution of ther- 
mal and quantum noise by examining the phase cor- 
relation function Cg{z — z') = (exp[i(6(z) — 6(z'))]) = 
exp[— A(\z — z'\)\. The function A that governs the de- 
cay of correlations, can be represented as a sum of the 
quantum (q) and thermal (th) parts: A = Aq + -4th- In 
the case of uncoupled quasi-condensates (J = 0) both 
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FIG. 4: (Color online). The universal (dimensionless) func- 
tions ICAth and ICAq as a function of distance (in units of 
healing length) for uncoupled quasicondensates. Thin solid 
line: contribution of the quantum noise; dashed line: its log- 
arithmic asymptotics. Thick lines: contributions of the ther- 
mal noise; the curves are labelled by the respective values of 
k B T/(gniD) 



have the same dependence on the Luttinger parameter, 
they are proportional to /C _1 . Consequently, ICA q and 
ICAth are universal functions, depending on k^T ' j '{griin) 
and \z — z'|/Ch only. In Fig. @]we plot these functions. 

Analyzing both contributions *4th,q to Cg we find: At 
small distances dAth/d\z — z'\ — ) 0; for distances \z — 
z '\ ^ Chff^iD / '(^bT) we recover the linear asymptotics 



Ath ~ kt\z — z'\. In contrast Aq is linear in \z — z'\ up 
to the healing length (h, for larger distances we obtain 
the asymptotics A q « ^ln[8|z — z'|/(7r(h)], i-e. Cg oc 
|z-z'|- 1//C - 

Consequently thermal fluctuations become dominant 
at \z — z'\ > (hgniD I {k^T) — tt/(kt1C)- These estima- 
tions of the relative contribution of the quantum noise 
are also valid for coupled systems with lj 3> 7t/(kt/C). 

To conclude, we have studied the fluctuation properties 
in samples of interacting quantum degenerate ID bosons. 
In contrast to previous work 0, Q, EH , our approach al- 
lows also to investigate also tunnel coupled, phase-locked 
ID systems, and provides a clear distinction between con- 
tributions of fundamental quantum noise and thermal 
excitations. In addition we show that on length scales 
\z\ £ 7t/(kt/C), where the fluctuations are dominated by 
thermal excitations, these systems can be described to a 
very good approximation by a simple semiclassical model 
based on the spatial evolution of the relative phase ac- 
cording to an Ornstein-Uhlenbeck stochastic process with 
Gaussian phase fluctuations. We expect that one can find 
similar semiclassical models for many other quantum de- 
generate systems at finite temperature, such as ID spinor 



systems [27 1. 
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